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Abstract 

In the paper, the temperature associated with a dynamical spherically symmetric black hole 
or with a cosmological horizon is investigated from the point of view of a point-like detector. 
First, we briefly review the Hamilton- Jacobi tunneling method for a generic dynamical spher- 
ically symmetric space-time, and present two applications of the tunneling method. Then, we 
apply a well-known relativistic quantum theoretical technique, namely the Unruh-DeWitt de- 
tector formalism for a conformally coupled scalar field in a generic FRW space-time. As an 
application, for the generic static black hole case and the FRW de Sitter case, making use of 
peculiar Kodama observer trajectories, the tunneling semiclassical results are fully recovered, 
automatically corrected by Tolman factors. Some remarks on the temperature of FRW universe 
are presented. For more general spaces interpolating de Sitter space with the Einstein-de Sitter 
universe a second set of poles is present, whose exact role remains to be clarified, plus an extra 
fluctuating term describing the way equilibrium is reached, similarly to de Sitter space. The sim- 
ple thermal interpretation found for de Sitter space is lost and forces, at a same time, a different 
quantum interpretation of the horizon surface gravity for the cosmological FRW models. 

PACS numbers: 04.70.-s, 04.70.Dy 



1 Introduction 

It is well known that Hawking radiation pQ is considered one of the most important prediction of 
Quantum Field Theory in curved space-time. Several derivations of this effect have been proposed 
[21 HI O [6] and recently the search for "experimental" verification making use of analogue models 
has been pursued by many investigators (see for example [3 E] ) . 
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In 2000 Parikh and Wilczek [9j(see also [10J) introduced the so-called tunneling approach for 
investigating Hawking radiation. Here, we shall firstly review a variant of their method, called 
Hamilton- Jacobi tunneling method [11 ^ I12| . [T3] and also |14j . This method is covariant and can be 
extended to the dynamical case [HI [161 EH 02] > an d to the study of decay of massive particles and 
particle creation by naked singularities [19] . In their approach, Parikh & Wilczek made a clever 
use of the Painleve stationary gauge for 4-dimensional Schwarzschild black hole 



ds 2 = -(l- ^\ dT 2 + + dr 2 + rW 



(1.1) 



which is regular on the trapping horizon rjj = 2mG 0. This is one of the key points since the use 
of singular gauges, as the Schwarzschild gauge, leads, in general, to ambiguities and it is useless in 
the dynamical case which we are interested in. 

The second merit we can address to the Parikh & Wilczek work was a treatment of the back- 
reaction on the metric, based on energy conservation. In the following, we shall limit to leading 
term results and neglect the back-reaction. However, it may be worth to recall that in the limit 
where the number of emitted quanta is reasonably large back-reaction effects can be accounted for 
by assuming that the mass parameter m is a continuous function of time T. Penrose's diagrams 
for this more general case have been determined too, e.g. in [20|, 121] , 

The Hamilton-Jacobi method is reasonably simple, even though subtleties are present (see, for 
example |22j). It is based on the computation of the classical action / along a trajectory starting 
slightly behind the trapping horizon but ending in the bulk, and the associated WKB approximation 
(c=l) 

Amplitude cx e 1 '* . (1.2) 

For an evaporating black hole such a trajectory would be classically forbidden since at the trapping 
horizon photons are only momentarily at rest, whence dr/dT < inside the horizoiH. The related 
semiclassical emission rate reads 

r x |Amplitude| 2 cx e~ 2 ^ . (1.3) 

with 9 standing for the imaginary part. In the tunneling across the horizon, the imaginary part of 
the classical action I stems from the interpretation of a formal horizon divergence, and in evaluating 
it, one has to make use of Feynman's prescription related to a simple pole in integration such as 

f dx^^- = »7r/(0) + real stuff. (1.4) 
x J _ x — tO 

This corresponds to the choice of suitable boundary conditions in quantum field theory approach. 

We may anticipate that in the WKB approximation of the tunneling probability, one asymp- 
totically gets a Boltzmann factor, in which an energy oj appears, i.e. 

Toce-K (1.5) 



lr The global event horizon is ro < ru for an evaporating black hole and ro > ru for an accreting black hole, as 
can be seen from the equation of radial null rays, ro = 1 — \/ru/ro- 
2 It means the photon must go back in time T to escape the horizon. 
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It is crucial in our approach that the argument of the exponent be a coordinate scalar (invari- 
ant quantity), since otherwise no physical meaning can be addressed to T. In particular, in the 
Boltzmann factor, f3 and the energy uj have to be separately coordinate scalars, otherwise again no 
invariant meaning could be given to the quantity j-. In the static case, we interpret T = h as the 
horizon temperature. 

In the cosmological case on the other hand, one can still have a trapping horizon despite the ab- 
sence of collapsed matter simply as a result of the expansion of the universe. As for an excreting 
black hole, this too is represented by a time-like hypersurface. Similarly, an approximate notion 
of temperature can be associated to such horizons based on the existence of a surface gravity and 
again the tunneling method gives a non vanishing amplitude having the Boltzmann form. However 
we will see that a comoving monopole detector seems to react to the expansion in a different, non 
"Boltzmannian" way, while reaching thermal equilibrium (or better, detailed balance conditions) 
only asymptotically for large times. 

The paper is organized in the following way: in Section 2, the Kodama-Hayward formalism and 
the Hamilton-Jacobi method are summarized; then, in Section 3, we consider the tunneling method 
for a generic static black hole in Kruskal-like gauge. In Section 4, the same formalism is applied 
to Friedmann-Robertson- Walker (FRW) spatially flat space-times. A formula for the de Sitter 
temperature which is valid in any coordinate patch is provided. Section 5, contains a discussion on 
the Unruh-DeWitt detector which is introduced in order to confirm through quantum field theoretic 
techniques the results of previous sections. In particular, a quite general formula for the response 
function per unit proper time is obtained and two relevant applications are presented. Section 6 
contains a further example, consisting of a (Q m ,£l\) cosmological model, where only approximate 
results can be obtained. Some toy models are presented and discussed. Concluding remarks follow 
at the end of the paper. 

We use the metric signature (— , +, +, +); Greek indices run over to 3 while mid-Latin as i, j 
only over and 1. We use Planck units in which c = H = G = kB = l- 



2 The Kodama-Hayward formalism and Hamilton-Jacobi tunnel- 
ing method 

In order to treat the spherically symmetric dynamical case, the use of invariant quantities plays a 
crucial role [17 \ \18 \ [23 | 124] . Here we review the general formalism. 

To begin with, let us recall that any spherically symmetric metric can locally be expressed in 
the form 

ds 2 = - ,,{.,■■ - 1 + R 2 (x i )dn 2 , i,j 6 {0, 1} , (2.1) 
where the two-dimensional metric 

dj 2 = lij {x i )dx l dx j (2.2) 

is referred to as the normal metric, {x 1 } are associated coordinates and R(x l ) is the areal radius, 
considered as a scalar field in the two-dimensional normal space. A relevant scalar quantity in the 
reduced normal space is 

X(x) = ^ (x)d i R{x)d ] R(x) , (2.3) 
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since the dynamical trapping horizon, if it exists, is located in correspondence of 

X (x) = 0, (2.4) 

provided that dix\n 7^ 0. The Misner-Sharp gravitational energy is defined by 

E MS (x) = ^R(x) [1 - X (x)] ■ (2.5) 

This is an invariant quantity on the normal space. Note also that, on the horizon, Ems\h = \Rh = 
M. Furthermore, one can introduce the Hay ward surface gravity associated with this dynamical 
horizon, given by the normal-space scalar 

1 

-[ 

2 



(2.6) 

H 



Recall that, in the spherical symmetric dynamical case, it is possible to introduce the Kodama 
vector field K, with (/C a G Q/3 ) ;/3 = that can be taken as its defining property. Given the metric 
(|2.ip . the Kodama vector components are 

V(x) = — L= e lj djR , K, 6 = = . (2.7) 

We may introduce the Kodama trajectories, and related Kodama observer, by means of integral 
lines of Kodama vector 



dx i i 1 
As a result 



V = ^=e ij djR. (2. 



dX dX y/— 7 

so, we have proved the following 

Lemma: In generic spherically symmetric space-times, the areal radius R is conserved along 
Kodama trajectories. 

In a generic spherically symmetric space-time a geometric dynamical identity holds true in 
general. This can be derived as follows. Let us introduce the normal space invariant 

T (2) = ^ijrp.. _ ( 2 1Q ) 

which we shall term reduced stress energy tensor trace. Then, making use of Einstein equations, it 
is possible to show that, on the dynamical horizon, (see, for example, 



K H = 777^ + ^R h t£ ] . (2.11) 
2R H 

Introducing the horizon area and the (formal) three-volume enclosed by the horizon, with their 
respective differentials 

Ah = ^R 2 H » dA H = 8irR H dR H , (2.12) 
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V H = -irR 3 H , dV H = 47rR 2 H dR H , (2.13) 



one gets 



^dA H = d^-f) + ^dV H . (2.14) 

This equation can be re-cast in the form of a geometrical identity, once we introduce the Misner- 
Sharp energy at the horizon [24J: 



kh j ( Ah \ T H 



(2) 



dM= t d {-f)-^t dV »- ^ 

We conclude summarizing the main ingredients of the Hamilton-Jacobi tunneling method. The 
Kodama vector, introduced above, gives a preferred flow of time and in this sense it generalizes 
the flow of time given by the Killing vector in the static case. As a consequence, we may introduce 
the invariant energy associated with a particle of mass m by means of the scalar quantity on the 
normal space, the Kodama, or generalized Killing energy, 

cu = -JC i d i I, (2.16) 

where I is the particle action which we assume to satisfy the reduced Hamilton-Jacobi equation 

f j d i Id j I + m 2 = 0. (2.17) 

As we allow for non-minimal gravitational coupling, the substitution m 2 — > m 2 + £JZ is in order 
whenever £ ^ 0, 1Z being the Ricci curvature scalar and £ a dimensionless coupling constant. 

If we are interested in Hawking effect, we may neglect the mass. One can reconstruct the action 
for particles coming out of the horizon by 



/ = dtdj+l drd r I (2.18) 

upon solving the Hamilton-Jacobi equation (|2,17p with zero mass. Then, one applies the near 
horizon approximation and the null horizon expansion which forces us to make use of regular 
gauges across the horizon. These assumptions allow one to know the classical action and compute 
its imaginary part making use of Feynman's prescription. The semiclassical result for the tunneling 
probability is then p3lfT7l [T8] 

r ~ e~ 23/ ~ e -Pn"H _ (2.19) 

In the static case, there is no doubt that Tjj = can be interpreted as the Hawking temperature; 
but what about the dynamical spherically symmetric case? From Einstein's equations and in 
presence of a dynamical horizon the geometrical dynamical Law holds as we have already seen. If 
also the Area Law for the black hole entropy is true, namely that 

S„-£, (2.20) 

then, we can rewrite the geometrical identity (|2,15p as a First Law of Thermodynamics for black 
holes, 

T (2) 

dhl = T H dS H - ~^-dV H , (2.21) 
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with 



T H = h^-. (2.22) 

Here we have re- introduced h in order to stress the quantum nature of entropy and temperature. 
Thus, it seems suggestive to interpret Th as a dynamical temperature of a slowly changing dy- 
namical black hole, interpretation which is also supported by the fact that the Hay ward dynamical 
surface gravity is an invariant for a generic spherically symmetric space-time. This interpretation 
has been put forward in [TTJ, [THJ Q21 [22]. Later in the paper, we shall try to present quantum 
theoretical arguments in order to substantiate this interpretation. 



3 Generic static black hole space-time 

As a first application of the formalism, let us consider a generic static black hole space-time. The 
starting point is a black hole metric in the Schwarzschild static gauge, 

dr 2 

ds 2 = -V(r)dt 2 + ^— + r 2 dtf , (3.1) 
W(r) 

where, for sake of simplicity, we shall assume W(r) = V(r), with V(r) having just simple poles in 

order to describe what we might call a nice black hole. Let ru be the (greatest) solution of V(r) = 0, 

the general formalism tells us that the horizon is located at r = r#; the Kodama vector coincides 

with the usual Killing vector (1,0,0,0); and the Hayward surface gravity is the Killing surface 

V' 

gravity, namely kh = K = ■ This gauge is singular on the horizon, and it is not appropriate for 
the Hamilton-Jacobi tunneling method. For this reason we now introduce the Kruskal-like gauge 
associated with this static black hole solution. The first step consists in introducing the tortoise 
coordinate 

Then one has — oo < r* < oo and 

ds 2 = -V{r)dt 2 + + r 2 dQ? = V{r*){-dt 2 + (dr*) 2 ) + r 2 (r*)dQ 2 . (3.3) 

Introduce Kruskal-like coordinates, according to 

R = - e Kr * cosh(/rf) , T = -e Kr * sinh(/rf) , (3.4) 

K K 

one has 

1 

K 



T 2 + R 2 = —e 2Kr , (3.5) 



and 



ds 2 = V(r*)e~ 2Kr *{-dT 2 + dR 2 ) + r 2 {T,R)dn 2 
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where now the coordinates are T and R, r* = r*(T,R), e*( r *) = V(r*)e~ 2Kr , and the normal 
metric turns out to be conformally flat. In this gauge, the metric is a spherically symmetric but 
time dependent one. The general formalism tells us that the horizon corresponds to 

(d T r) H = ±(d R r) Hl (3.7) 

and this is equivalent to T = ±R and r* — > — oo. The Killing-Kodama vector is 

JC = e -*^ {d R r,-d T r) . (3.8) 

Making use of the general formula, a direct, but tedious calculation confirms that the Hayward's 
surface gravity is still the Killing one. In fact, we have 

KH= e ^{-d 2 T r + dlr) H = ^-. (3.9) 

Let us apply the Hamilton-Jacobi tunneling method, working in this time-dependent gauge. 
The Kodama energy is 

w = e - * (r * ) (d R r&rI-d T rd R I) (3.10) 

where I is the classical action. In this diagonal conformally flat gauge, the Hamilton-Jacobi equation 
is simply dxl = ±d R I and the null horizon expansion condition gives dT = ±dR. Making the + 
choice (outgoing particle), in the horizon approximation, one has 



/ ~ 2 J dRd R I. (3.11) 
Eq. ([3TTUD gives 

9rI = a e \ • (3.12) 
o R r — ot r 

The near horizon expansion gives 

d R r - d T r ~ (8 2 R r - d$r) H (R - R H ) . (3.13) 
so that the action (|3.1ip becomes, 

dR—- , (3.14) 

K H (R-R H -iti) v ; 

Making use of Feynman's ie— prescription for the simple pole, we finally get 

S/=^ = ^. (3.15) 

This imaginary part can be interpreted as arising because of a non-vanishing tunneling probability 
rate of (massless) particles across the event horizon, 

r ~ exp(-23f/) ~ e~v^> WH . (3.16) 

V 

The well-known result Tjj = is recovered. 
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4 Cosmological horizons 



As a second application of the formalism, we consider a generic FRW space-time with constant 
curvature spatial sections. Its line element can be written as 



ds 2 



-dt 2 +a 2 {t) 



dr 



1 — kr 2 



+ [a(t)r] 2 dn 2 . 



(4.1) 



Here k := p, where I is such that a(t)l is the curvature radius of the constant curvature spatial 
sections at time t and, as usual, k = 0,-1, +1 labels flat, open and closed three-geometries, 
respectively. In this gauge, the normal reduced metric is diagonal and 



X (t,r) = 1- [a(t)r]' 



H 2 (t) + 



a 2 {t) 



The dynamical trapping horizon is implicitly given by xh = 0, namely 



R H ■= a(t)r H 



1 



with 



H(t) 



(0 



d(t) 
a(t) ' 



(4.2) 



(4.3) 



provided the space-time energy density p(t) is positive. It coincides with the Hubble radius as 
defined by astronomers for vanishing curvature, but we shall call it Hubble radius in any case. The 
dynamical surface gravity is given by equation (12. 6p and reads 



H 2 {t) + ~H(t) + 



k 



2a 2 (t) 



R H (t), 



(4.4) 



and the minus sign refers to the fact the Hubble horizon is, in Hayward's terminology, of the inner 
type. According to (|2.7p . the Kodama vector is 



JC = \/1- kr 2 (d t -rH(t)d r ) 
so that the invariant Kodama energy of a particle is equal to 



UJ 



1 — kr 2 (—dtl + rH{t)d r I) = \J I - kr 2 u 



(4.5) 



(4.6) 



Notice that K. is space- like for ra > (H 2 + fc/a 2 ) -1 / 2 , i.e. beyond the horizon. It follows that we 
can only ask for particles to be emitted in the inner region r < rn- 

The next ingredient is the reduced Hamilton-Jacobi equation for a relativistic particle with mass 
parameter m, 



(d t I) 2 + { -^l(d r I) 2 +m 2 



0. 



(4.7) 



Making use of (|4,6p . one can solve for d r I, namely 



aHCo(ar) ±a ( w 2 



d r I 



m 



2 + m 2 ( H 2 + A 



(ar) ; 



H 2 + 



(ar) 2 



(4.8) 
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with the signs chosen according to which direction we think the particle is propagating. The 
effective mass here defines two important and complementary energy scales: if one is interested in 
the horizon tunneling then only the pole matters (since the denominator vanishes), and we may 
neglect to all the extents the mass parameter setting m = (since its coefficient vanishes on the 
horizon) . 

On the opposite, in investigating other effects in the bulk away from the horizon, such as the decay 
rate of composite particles, the role of the effective mass becomes relevant as the energy of the 
particle can be smaller than the energy scale settled by m, and the square root can possibly acquire 
a branch cut singularity. 



4.1 The FRW space-time 

As an application of the last formula we may derive, following [18] (see also [25J), the cosmic horizon 
tunneling rate. To this aim, as we have anticipated, the energy scale is such that near the horizon, 
we may neglect the particle's mass, and note that radially moving massless particles follow a null 
direction. Then along a null radial direction from the horizon to the inner region we have 

At = a{ t\ Ar. (4.9) 

V 1 — kr 2 

The outgoing particle action, that is the action for particles coming out of the horizon towards the 
inner region, is then 



/ = J dtd t I + J drd r I (4.10) 
= 2 / drd r I (4.11) 



upon solving the Hamilton- Jacobi equation (|4.7p with zero mass and using (|4.9p . For d r I we use 
now Eq. (14. 8h . which exhibits a pole at the vanishing of the function F(r,t) := 1 — (a 2 H 2 + k)r 2 , 
defining the horizon position. Expanding F(r, t) again along a null direction, one gets 

F{r, t) « +4n H a{t){r - r H ) + . . . , (4.12) 

where kh is the Hayward dynamical surface gravity. In order to deal with the simple pole in the 
integrand, we implement Feynman's ie - prescription. In the final result, beside a real (irrelevant) 
contribution, we obtain the following imaginary part |18] 



31 = -^. (4.13) 

This imaginary part is again interpreted as arising because of a non- vanishing tunneling probability 
rate of (massless) particles across the cosmological horizon, 

r ~ exp(-23I) ~ e~*%5J UB . (4.14) 

Notice that, since kjj < and loh > for physical particles, (|4.13p is positive definite. As showed 
in [18], this result is invariant since the quantities appearing in the imaginary part are manifestly 
invariant. Furthermore T = — k#/27t satisfies a First Law. As a consequence, at least, in some 
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asymptotic regime and for slowly changes in the geometry, we may interpret T = —kh/^tt as the 
dynamical temperature associated with FRW space-times. 

In particular, this gives naturally a positive temperature for de Sitter space-time, a long debated 
question years ago, usually resolved by changing the sign of the horizon's energy. It should be noted 
that in literature, the dynamical temperature is usually given in the form T = M- (exceptions are 
the papers [261 EZ]) with H 2 = A/3 = Hq. Of course this is the expected result for de Sitter space 
in inflationary coordinates, but it ceases to be correct in any other coordinate system since, for 
example, H = Hq tanh Hot in global coordinate system with positive spatial curvature. In this 
regard, the H and k terms are crucial in order to get an invariant temperature. Since this fact 
seems not so widely known, for sake of completeness, we shall try to show it in greater detail. 

de Sitter space in the global patch is described by the metric 



ds 2 = -dt 2 + cosh 2 (H t) 



(4.15) 



.1 — kr 2 

with a(t) = cosh(-ffoi)) and k = Hq. The Hubble parameter is time dependent 

H(t) = H tanh{H t) , (4.16) 
and satisfies the identity 

H(t) = H 2 -H 2 (t) = J^. (4.17) 

Making use of (|4.3|) . the horizon radius is 

R H := a(t)r H = . 1 = tt (4.18) 

v^w + A) ° 

as it should be, and the Hayward's surface gravity 

K H = ~ (^H 2 (t) + iff(t) + R H (t) = H , (4.19) 

as it should, since it is an invariant quantity. Hence, we see that the H and k terms have to be 
present in a generic FRW space-time. The important spatially flat case straightforwardly follows, 

Note that this is independent on position, suggesting that kh really is an intrinsic property of FRW 
space linked to the bulk. The horizon's temperature and the ensuing heating of matter was foreseen 
several years ago in the interesting paper [28]. We are going to challenge this interpretation. 
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5 Quantum thermometers 



We recall that for the decay rate of a massive particle in de Sitter space, the exact quantum 
field theory calculation of Moschella et al. [29] supports the WKB semiclassical tunneling result of 



What about the other energy scale, the one associated with the horizon tunneling? We have 
shown that the semiclassical WKB method leads to an asymptotic particle production rate, involv- 
ing a "temperature" 



where kh is the Hayward invariant surface gravity. For a generic spherically symmetric space-time, 
this result obtained by the Hamilton-Jacobi tunneling method seems a very clear prediction, namely 
an answer to the question: how hot is our expanding universe? A possible way to understand this 
issue using quantum field theory in curved space-time is to make use of a "quantum thermometer" 
(basically, an Unruh-DeWitt detector) and evaluate its response function, that is, loosely speaking, 
the number of clicks per unit proper time it detects as it is carried around the universe. For a 
recent review, see [31j . 

In our approach, since we would like to obtain an invariant result, we will consider detectors 
which follow Kodama trajectories in a generic spherically symmetric space-time. The problem of 
back-reaction on these Kodama trajectories has been investigated in |32j. 

As we will see, the Unruh-DeWitt thermometer gives a clean answer only in the stationary 
case, and for FRW space-time the situation is not so simple, since horizon effects are entangled 
with highly non trivial kinematic effects. For general trajectories in flat space-time see the recent 
paper [33] , 

An interesting analysis has been also put forward by Obadia [34] . In a recent paper [35], local 
scaling limit techniques have been used in investigating the Hawking radiation. 

5.1 Unruh-DeWitt detector 

We recall that the Wightman function of a free, conformally coupled, massless scalar field (the only 
case we shall be interested in here) may be written as 



[301 033 • 







a 



where the modes functions f a ( x ) satisfy the equation 




(5.3) 



As it is well known, it is convenient to introduce the conformal time r\ defined by 



(5.4) 



Thus, the flat FRW space-time becomes conformally flat, 




(5.5) 
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Let us denote by x = (t],x). Making the ansatz 

/ s (x) = ^e-**, (5.6) 

one has, for the unknown quantity g(k), 

g"(v) + k 2 g(ri) = . (5.7) 

As a consequence, making the choice of the vacuum given by 
ir)\k\ 



one formally has 

W(x,x') = / 4 -— — - . (5.9) 

This expression is meaningless as it stands and it has to be interpreted in the sense of distributions. 
Within Unruh-DeWitt detector issue, the usual prescription reads, [3], 

/dk piiv—v' — ie )\k\—"ik-(x-x') 
—r , , / / \ • (5.10) 

2k a{rj)a{rj') V ; 

Integrating on k, one arrives at 

W{x,x') = lim — — — — . (5.11) 

e-s-0 47r z ayq)a{rf) \x — x \ — \r] — r( — ley 

However, it has been shown by Takagi [36] and Schlicht [37] that this prescription is non manifestly 
covariant. Since one is dealing with distributions, the limit e — > + has to be taken in the weak 
sense, otherwise it may lead unphysical results with regard to instantaneous proper-time rate in 
Minkowski space-time for a Unruh-DeWitt detector. 

Here, we adapt the Schlicht proposal (see also |38} 139]) to our conformally flat case, namely 

/rlk p—ik(x—x')—ie(x+x') 
—r --— , (5.12) 

2k a(rj)a(rf) ' v 1 

where an over dot is for derivative with respect to proper time. Integrating on k, one arrives at 

W(x,x')= lim — ^ - r T . — — • (5.13) 

e->0+ iir z a{rj)a{rj'j [{x — x j + ie{x + x )\ z 

In the flat case, this result has been generalized by Milgrom and Obadia, who made use of an analyt- 
ical proper-time regularization [39] . It should be noted the appearance of Minkoswki contribution, 
as a function of the conformal time rj. 

The transition probability per unit proper time of the detector depends on the response function 
per unit proper time which for radial trajectories may be written as 

dF 1 f T ~ T ° e -iEs 

— = lim — -Re / ds rr~r^ 7 s T~r\ 7 ^TTo > (5-14) 

dr e-s>o+ 2tt z Jq a{T)a{T — s)[x(t) — x(t — s) + ie(x{T) + x(t — s))\ z 
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where tq is the initial time, and E is the energy associated with the excited detector state. In 
some cases one can take the limit To — > — oo, but not in general since as a rule there is a Big Bang 
singularity. Instead we will be interested in the long term behavior of the detector, that is in the 
large r limit at fixed tq. The covariant ie-prescription is necessary in order to deal with the second 
order pole at s = 0. One may try to avoid the awkward limit e — > + by omitting the e-terms but 
subtracting the leading pole at s = (see for details [10]). In the calculation, the normalization 
condition 

g^x v = [a(r)x(r)] 2 = -1 (5.15) 

characteristic of time-like four-velocities plays a crucial role. As a consequence, for At > 0, after 
some calculations, one can present the detector transition probability per unit time in the fornJl 



dF_ 

~dr ~ 2vr 2 



1 r° , x ( 1 1 \ 1 f°° , cos(Es) 

- 2 j o dsco S (Es)^-^ )+ - 2 j-— 2 J^ds- T{ - j , (5.16) 



where now there are no ie-terms and 

o- 2 (r, s) = a(r)a(r - s)[x(t) - x(t - s)] 2 . (5.17) 

It should be noted the presence of the last fluctuating finite time tail term. In some cases it controls 
how fast the thermal equilibrium is reached where, to recall the reader, At is the duration of the 
experiment, i.e. the time the detector is switched on. 

We can go further, observing that due to (|5.15p . for small s, one has 

a 2 (T,s) = -s 2 [l + 0(s 2 )}. (5.18) 

Thus we may write 

ct 2 (t,s) = -s 2 [l + s 2 d(T,s)] . (5.19) 

As a result, we may introduce the even <t 2 (t, s) and odd part <j 2 (t, s) in s of the quantity ct 2 (t, s) 
and arrive at the two contributions 

dF dF e dF Q 1 f°° , cos(Es) 

" + (5 - 20) 

where 

dF e 1 f°° , , s ( 1 1 \ , 

■d7 = 4^L dscosiEs) {^) + 7*)' (5 ' 21> 

« > f 4 ?S!S. (5.22) 

This is the main formula which we will use in the following. One can view the pole subtraction as 
an elementary example of a renormalization procedure, one that would not be necessary where the 
e-terms were kept all along. However, Eq.s (|5.16p . (|5.20p are much more convenient to deal with 
than the original expression containing the e-terms, since the limit in distributional sense must also 
be taken at the end of any computation. 



This is not quite the original expression found in [40] but can be obtained from it by simple manipulations. 
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In the important stationary case in which a(r,s) 2 = <J 2 (s) = a 2 (—s), for E > and At > , the 
odd part drops out and one simply has 

iF - 1 r^w i + n_ i r*^. „ 23) 



dr Air 2 J_ 00 \<t 2 (s) s 2 J 2tt 2 J^ t <j 2 (s) 

In these cases, examples are the FRW de Sitter space and static black holes, we shall take the limit 
Ar —> oo, and the fluctuating tail vanishes. 

5.2 The generic static black hole revisited 

It is clear that the Unruh-DeWitt detector formalism can be applied to the static black hole of 
Section 3. The key point to recall here is that in the Kruskal gauge (j3.6|) the normal metric 
is conformally related to two dimensional Minkoswki space-time, and the normal metric is the 
important one for radial trajectories. The second observation is that the the Kodama observers are 
defined by the integral curves associated with the Kodama vector, thus the areal radius r(T,R) is 
constant, and r* is constant. As a consequence, the proper time along Kodama trajectories reads 

dr 2 = V{r*)dt 2 = e^^ ] {dT 2 - dR 2 ) = a 2 {r*){dT 2 - dR 2 ) , (5.24) 



so that 



t = , , (5.25) 



and 



R( T ) = - e Kr * cosh ( k - T ) , T(r) = -e Kr * sinh ( k - T ) . (5.26) 

Recalling that the geodesic distance is 

a(r, s) = e*^) [- (T(r) - T(r - s)) 2 + (R(r) - R(r - s)) 2 ] , (5.27) 

one gets for T(r) and R(t) of the form ([5.26P 

W(r*) -/ ks 
smh 2 — 
k 2 \2y/V(r* 

Since <t 2 (t, s) = cr 2 (s) = a 2 (—s), we can use the formula (I5.23[) in the limit r goes to infinity: 

dF k r°° , ( 1 1 \ 

dxe 1 « — —2— + -o , (5.29) 



(5.30) 



dn~ 8ir 2 VV* J-oo V sinh 2 x a; 2 . 

The integral can be evaluated by the Theorem of Residues and the final result is 
dF 1 E 



dr 2vr Wv^e 

e « — 1 



Since the transition rate exhibits the characteristic Planck distribution, it means that the Unruh- 
DeWitt thermometer in the generic spherically symmetric black hole space-time detects a quantum 
system in thermal equilibrium at the local temperature 

T = — (5.31) 
2nVV* 
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With regard to the the factor v V* = \f— goo, recall the Tolman's Theorem, which states that for 
gravitational system at thermal equilibrium, one has T^—goo = Constant. For asymptotically flat 
black hole space-times, one obtains the "intrinsic" constant temperature of the Hawking effect, i.e. 

= £ = < 5 - 32 > 

It is not necessary to stress how this result agrees with several derivations of Hawking effect, 
including the tunneling approach we have already discussed. 

We would like to point out that this is a quite general result, valid for a large class of nice black 
holes, as for example Reissner-Nordstrom and Schwarzschild-AdS black holes. On the other hand, 
the Schwarzschild-dS black hole cannot be included, due to the presence of two horizons. However, 
as an important particular case, we may consider the static de Sitter black hole, defined by 

V(r) = l-H 2 r 2 , #o = |- (5-33) 

The unique horizon is located at rjj = and Gibbons-Hawking temperature is [41J 

7i-§. (5.34) 

We will present another derivation of this well known result in another gauge in the next Section. 

We conclude this Section, making some remarks on de Sitter and Anti-de Sitter black holes. 
First, we observe that in a static space-time, namely the one corresponding to a nice black hole, 
the Killing-Kodama observers with r = K constant, have an invariant acceleration 

A 2 = g» u A»A» = , (5.35) 

where A^ = u^V^u^, being the observers 4- velocity. In the case of de Sitter black hole, one has 

A= -r^- <«■*> 

As a result, 



and the de Sitter local temperature felt by the Unruh detector, 
H Q 1 



'-dS 



2tt y/1 - HqK 2 

can be re-written as 
1 



(5.38) 



Tds = ^ ^A 2 + H 2 = ^T 2 + T 2 H . (5.39) 

A similar result was also obtained for AdS in |43j . and it reads 
1 



T AdS = ^A 2 -H 2 . (5.40) 
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We would like to show that it is again a particular case of our general formula (j5.31|) . In fact, 
it is sufficient to apply it to the four-dimensional topological black hole with hyperbolic horizon 
manifold found in |44l li"5l l4"6l |4"T] , which is a nice black hole with 

C 

V(r) = -1 + H^r 2 , (5.41) 

where C is a constant of integration related to the mass. The space-time is a solution of Einstein 
equation with negative cosmological constant A = —3Hq, which is asymptotically Anti-de Sitter. 
When the constant of integration goes to zero, one has still a black hole solution, and calculation 
similar to the one valid for de Sitter space-time gives 



Hp 1 = I 

2tt yZ-l + H'^K 2 2tt 



T AdS = ^ - . | =^= = - X /A2 - Hi , (5.42) 



which is the Deser et al. result. Thus, for symmetric space time with constant curvature one has 
that the local temperature felt by the Kodama-Unruh-DeWitt detector can be written as 



T=y/1% + ,T 2 H , (5.43) 

where Tjj is the Unruh temperature associated with the acceleration of the Kodama observer, Tqh 
is the Gibbons-Hawking temperature and ? = 1 for the de Sitter space-time, ? = 0, for Minkowski 
space-time (this is the original Unruh effect) and q = — 1 for the "massless" AdS topological black 
hole. This formula may help to understand better the relation between the Unruh-like effects and 
the genuine presence of a thermal bath and shows that, in general, the Kodama-Unruh detector 
gives an intricate relation between Killing-Hayward temperature and other invariant temperatures 
such as the Unruh's one. 



6 The FRW and de Sitter expanding universes 

To apply the Unruh-DeWitt detector formalism to cosmology we consider a generic FRW spatially 
flat space-time. This case has been investigated also in [34J. Recall that here the areal radius is 
R = ra(t). Thus, for the Kodama observer, one has 

r(t) = 4v (6- 1 ) 
a(t) 

with constant K. For a radial trajectory, the proper time in FRW is 

dr 2 = a(i])(dri 2 - dr 2 ) . (6.2) 
As a function of the proper time, the conformal time along a Kodama trajectory is 

ti(t) = — I dr dri 

J a{ri)y/l - K 2 H 2 (t) 

dr , (6.3) 
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H(t) being the Hubble parameter as a function of proper time. In general, we may use Equation 
(|5.20p in which, for radial Kodama observer, one has 

x(r) = (r ? (r),r(r), 0,0)= f- / \ dr,^- 0,o] . (6.4) 

\ J a[r)y/V{T) a(r) J 

A very important example of FRW space is the stationary flat de Sitter expanding (contracting) 
space-time, which, in the FRW contest, is defined by a(t) = e Hot . Thus, 

ds 2 = -dt 2 + e 2Hot {dr 2 + r 2 dQ 2 ) . (6.5) 

Here H(t) = Hq is constant as well as V = Vq = 1 — H^K 2 . For Kodama observers 

r= V / V r t, a(r)=ev^ , (6.6) 

and 

l -4a= T -4a= T 

r](T) = -—e V% , r(T) = Ke V^o , (6.7) 
Ho 

so, the geodesic distance is 



^ = --^nh^j . (6.8) 

This result is formally equal to the one obtained for the generic static black hole (|5.28p . Since again 
<t 2 (t, s) = cr 2 (s) = a 2 (—s), we may use (|5.23p and obtain, for E > and in the infinite time limit 

— = " = / dxe H o — + — . 6.9) 

dr 8V%k 2 Loo V sinh 2 x x 2 ) V ; 

Again, we arrive at 

dF 1 E 



dr 2tt 2 ^V^q e 

e H o - 1 



3.10) 



which shows that the Unruh-DeWitt thermometer in the FRW de Sitter space detects a quantum 
system in thermal equilibrium at a temperature T = 2 ^jy Here, the Tolman factor is substituted 



by a Lorentz factor, which represents, as we already know, the Unruh acceleration part. In fact we 
recall that 4-acceleration of a Kodama observer in a FRW space-time turns out to be 



A 1 = A^Afj, = K 



H{t) + {l-H 2 (t)K 2 )H 2 {t) 



(1 - H 2 {t)K 2 ) 



(6.11) 



where A^ := u^V^u^, being the 4-velocity of the detector, that is the (normalized) tangent 
vector to the integral curves of the Kodama vector field. As a result, for dS space in a time 
dependent patch we have 

,2 K 2 Hf: 
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showing that 



7/11 'm+A* 



in agreement with the dS static calculation. When K = 0, one has Vq = 1 and the classical 
Gibbons-Hawking result T^s = ^ is recovered. 

We close this section with a brief discussion of an asymptotic phenomenon which will be relevant 
to the following sub-section, first noted in de Sitter space by Garbrecht, Prokopec [48]: how is the 
thermal distribution of the response function reached in the limit of very large times? To this aim 
we consider the finite time contribution due to the fluctuating tail (the last term in Eq. (15.23P ) for 
de Sitter or Black hole cases. A direct calculation of the tail using Eq. (|5.28p for black holes or 
Eq. (USD for dS, gives 



K 2 f°° , COs(Es) K ™ ne -2™TflAr / £ \ 

L d '^f = *?£"^3T fe- cos(EAT) - sm(EAr) ) (6 - l3) 

where k is the horizon surface gravity and Tjj is the local Hawking temperature of the horizon, 
given by k/2-k^/V . We recall that k = Hq for de Sitter space and V (r#)/2 for the black hole. 
As a result, the fluctuating tail term drops out exponentially. For large E/Th the oscillations 
are at first very large relative to the equilibrium value and for few Hubble times then they decay 
moderately fast, reaching values comparable to equilibrium values only at times At ~ a few E/Th- 
For instance with E/Th = 15 we have a plo1@ 




although only at At ~ 90 is the tail smaller than the equilibrium value, which if of the order 
10 -41 . From this point of view the oscillating terms seem more sensible than the equilibrium value. 



6.1 Asymptotically de Sitter space-times 

We end with some considerations about the non stationary case which appear when space-time is 
only spatially homogeneous as in generic FRWL models, but limiting our computations to comoving 

4 Obtained using Wolfram's Mathematica numerical code. 
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detectors (i.e. K = 0), the less interesting from the point of view of the Unruh effect. It should 
be clear from previous considerations, that the formula for the transition rate of the detector 
(|5,16p is manageable, in the sense that we can extract quantitative information, only in few highly 
symmetrical circumstances. For instance the case of a realistic cosmology dominated by matter 
and vacuum energy has a scale factor of the form 

/o \2/3 / q \ 1/3 

a{t) = a smhl-^hiH t\ , a =(^\ (6.14) 

where Hq is the present value of the Hubble constant. We put h = 3\/QxHq/2 for simplicity: then 
the conformal time is given up to a constant by the hyper geometric function 

— ^ = sech 2/3 (M) 2 Fi(5/6,l/3;4/3;sech 2 (/it)) + C (6.15) 

a[t) a h 

where C is possibly a complex constant but the first term is real. Clearly integrals such as (|5,16p . 
(|5.2ip and (|5.22p are now outside an analytical treatment and one can resort to numerical eval- 
uations. We only note few things: at small ht -C 1 the model follows the Einstein-de Sitter law 
a(t) ~ t 2 / 3 of a matter dominated universe. In this case the geodesic distance is 



<t 2 (t, s) = 6t(t — s) sinh 2 



St 1 "" 
6 V r 



and the pole structure is the following: s = r with residue (also a branch point) and s = 0, 
a double pole with residue —2iE. So this detector will not take notice of the expansion of the 
universe, even if the horizon temperature is non vanishing (T# ~ t~ 1 )- Second, there will be in 
general further poles over those along the imaginary axis, due to the factor a(r — s) which multiplies 
the geodesic distance and which may have zeroes when extended throughout the complex plane: 
for example with a(r) as given we have further poles at 

mri 

Sn=T + —— 

h 

which however seem to disappear in a large time limit. Moreover, we have some numerical evidence 
that their effect is eaten by the oscillating tail that we discussed in the context of dS space. 
Finally, we have the noteworthy feature that the geodesic distance appearing as integrand in the 
response function is time dependent, which means the detector will not feel a thermal bath except 
asymptotically at large Hubble times. In fact at large ht 3> 1 the scale factor describes de Sitter 
space with the smaller Hubble constant \/Q\Hq and <t 2 (t, s) exhibits the characteristic dS double 
poles at s n = 2mri/^0,\Ho, so we expect an asymptotic thermal distribution with T = ^JVL\Hq/2tt. 
This can be seen in few simple toy models, which are not physically motivated, but that allow some 
analytical evaluation. For instance one could try the scale factor a{t) (always with spatial curvature 
k = 0) 

a(t) = cosh 2 (^ft) ■ ( 6 - 16 ) 

This scale factor tends to a de Sitter expansion phase at late times and it is singularity free. The 
Unruh-DeWitt detector is understood to sit at the origin of the comoving mesh so we have the 
huge simplification 

r(i) = t. (6.17) 
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a\r,s) = -o(r)o(r - s) [ V (r) - r,(r - s)Y = --^sinh^ I -^t ) = <r^( S ) . (6.19) 



As a consequence, the conformal time is 

T/(*) = ~tanhf^ , (6.18) 
and the geodesic distance reads 

* sinh* (^ - -» 

Clearly the transition probability rate is given by the equation (|5.20p . for E > 0, and it is formally 
equal to the de Sitter case. We stress that this is only true for Kodama trajectories, r(t) = K/a(t), 
with K = 0. For non vanishing K, ij(t) may be expressed as an intractable elliptic integral, and 
the transition probability rate is time dependent. The same result can be obtained by the choice 

a(t) = sinh 2 (Jy*) ' ( 6,2 °) 

this toy model being asymptotically de Sitter, but close to Milne model and singular near t = 0. 

The following toy models are less trivial but, again for K = 0, they exhibit features which are 
generally present for more realistic models. They are described by the choices 

cosh(f/oi) 

and 

a ^ = ■ urr A ■ ( 6 - 22 ) 

These models represent asymptotically contracting dS spaces, and their geodesic distance reads 
respectively 

2 / \ 2 / \ cosh 2 H (t - s/2) 

a (r, s) = a dS (s) — — — — , (6.23) 

cosh Hot cosh Ho(t — s) 

2 f \ 2/\ sinh 2 H (t- s/2) 

a (t, s) = cr,o(s) . rT . - TT . r . (6.24) 

K J dby ' sinh H t sinh H (t - s) v ' 

We will investigate only the first one. Using trigonometric identities one can write 

1 _ 1 H 2 



(6.25) 



a 2 (t,s) a 2 s (s) 4 cosh 2 (H (t - s/2)) 
Thus the response function is the de Sitter one, which we already discussed, plus the integral 
Hi f At cos Es 

J = -tA ot ; ~-ds 6.26 

8vr 2 y cosh 2 (H (t- s/2) V ; 

Using 

oo 

sech 2 (x) = -4^(-l)\e 



2nx 

->-) " c 
n=l 
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a term by term integration gives 




n e 



-nH At - 



n 



n(cosEAt - 1) + — sin E At 
Ho 



(6.27) 



n 



2 + 




The similarity with de Sitter case, Eq. (|6.13p . is quite evident so some of the features of the pure 
dS case are indeed present. To summarize, we would say that the detector clicks close to a de 
Sitter response and reaches thermalization through decaying oscillations as r goes to infinity. In 
fact, as far as the regime Hqt 3> 1 is concerned, expression (|6.27p simplifies to the de Sitter space 
result. We may think of this as describing a de Sitter thermal noise continuously perturbed by the 
expansion (or contraction) of the universe. In particular, insofar as we can speak of temperature, it 
registers the de Sitter temperature, different from the horizon temperature parameter which in the 
present case is 2-kTh = Hq — 2Hq sech(2i7ot). This illustrates the feature that comoving detectors 
are probably unable to reveal the horizon temperature, whose quantum interpretation as an Unruh 
effect requires non trivial Kodama trajectories. But for such accelerated detectors {K ^ 0) one 
expects even more complicated results as Eq. (|5.43[) may lead one to think. 

Conclusions 

In this paper, making use of the Hamilton- Jacobi tunneling method, we have re-derived the Hawking 
effect for a generic static black hole and the analogous tunneling phenomenon which appear to exist 
also in a generic FRW space-time possessing a trapping horizon. The method also works for slowly 
changing black holes, except that the event horizon gets replaced by a time- like surface known 
as the Hayward's trapping horizon. With the aim to better understand the temperature-versus- 
surface gravity paradigm, the asymptotic results obtained by this semiclassical method have been 
tested with more reliable quantum field theory techniques as the Unruh-DeWitt detector analysis. 
For black holes the two analysis are mutually consistent and even predict the dependence of the 
temperature on position or acceleration. For cosmology and away from de Sitter space the thermal 
interpretation is lost but the detector still gets excited by the expansion of the universe: the second 
set of poles, as well as the presence of an odd part, arising in a quasi-de Sitter universe signals this 
lack of thermal behavior. Note, this result seems to show how the thermal interpretation breaks 
down in most of the cases: the time-dependence of the transition rate is expected to be persistent 
in all purely non-de Sitter dynamical solutions. Moreover, from the comoving detector point of 
view the horizon surface gravity seems less prominent and not unambiguously associated with a 
temperature parameter as the tunneling method suggests. For instance, in the Einstein-de Sitter 
regime there seems to be no excitations (for co- moving detectors only). It remains to see whether 
there is any non trivial effect on accelerated Kodama, or more general, trajectories. 
The fact that our co-moving monopole detector apparently seems unaware of the trapping horizon 
is related to the ambiguity of the particle concept in cosmology and forces, in our view, a different 
interpretation of the tunneling picture, hopefully giving matter to future work. One possibility is 
that the horizon surface gravity could represent an intrinsic property of the horizon itself, leading 
to some kind of holographic description, while the detector in the bulk simply clicks because it 
is immersed in a changing geometry. In fact, we would expect the clicks in almost any changing 
geometry, even for those lacking a trapping horizon. 
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